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Abstract: The survival probability for exclusive diffractive Higgs production is calculated. The contribu- 
tion of short distance interactions are taken into account, by summing over Pomeron loops in perturbative 
QCD. The summation is performed by developing an iterative technique to sum over loop diagrams with 
higher and higher generations of loops. The results show that the survival probability depends inversely 



L* on energy and is small for the LHC range of energies, and could be even less than 1 %. 

• rH 

X 



S— i Keywords: BFKL Pomeron, Triple Pomeron vertex, Higgs boson, summing Pomeron loops, QCD 



*Email: jeremymi@post.tau.ac.il; jeremy.miller@ist.utl.pt 



The most important event anticipated at the LHC is the detection of the Higgs boson, in proton 
proton collisions. The most desirable result is the exclusive production of the Higgs with no other particles 
produced in the scattering. As such there are large rapidity gaps (LRG) between the Higgs and the 
emerging protons shown in Fig. |l]. Thanks to the large rapidity gaps, exclusive Higgs production has 
an excellent experimental signature, and offers the best chance for successfully isolating the Higgs boson. 
Fig. [l] is a double t-channel gluon exchange leading to a zero net color flow in the t-channel, which 
cancels the possibility of additional inelastic scattering apart from the Higgs. The evolution of ladder 
gluons between the two t channel gluons forms the structure of the BFKL Pomeron in the leading log 
approximation. Unfortunately in high energy scattering, the production of extra unwanted particles from 
more parton showers, could spoil the large rapidity gaps as shown in Fig. g which means detecting the 
Higgs is problematic. The survival probability is the probability of having large rapidity gaps between the 
Higgs and the emerging protons, and as such its value characterizes the probability for exclusive Higgs 
production without further unwanted production. 




Figure 1: Exclusive Higgs production from t- Figure 2: The production of Higgs with extra pro- 
channel BFKL Pomeron exchange, with large rapid- duction which spoils the LRGs, arising from addi- 
ity gaps (LRG) between the Higgs and the emerging tional inelastic scattering. . 
protons. 

The contribution of the full set of parton showers that destroy the LRGs is required for estimates of 
the survival probability. Both long and short distance interactions contribute, so a reliable estimate is 
difficult to obtain. The contribution of short distance interactions can be treated in perturbative QCD, 
and for high density protons this type of hard process can be treated by summing over ladder diagrams of 
the type shown in Fig. |l[ In this approach Fig. p] is the sum over all diagrams with n rungs of the ladder. 
The vertical lines of the ladder are themselves a superposition of the sum over n rung ladder diagrams and 
so on. This sum over ladder diagrams leads to the scattering amplitude in the leading log approximation 
(LLA) which is proportional to [1-5] A (s, t) oc Yl^Lo ( a (*) m ( s / s o)) n /n\ = s a<yt \ where a (t) is the Regge 
trajectory. Hence the sum over ladder diagrams is achieved by replacing the pair of interacting vertical 
gluons in Fig. |l] by a "reggeon" which behave as s a ^ for energy s> 1. According to the optical theorem 
the total cross section behaves as a tot oc s a{ ^~ x . Experimentally it is known that the total cross section 
rises slowly with s which means a (0) > 1. Pomeranchuk first commented [6, 7] that this behavior is 
matched by the theoretical prediction that a (0) > 1 when the t-channel exchange carries zero quantum 
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numbers, including zero charge and color flow in the t-channel. Such particles with quantum numbers 
of the vacuum exist in QCD for bound gluon states. This trajectory is called the Pomeron named after 
Pomeranchuk, which is the double t-channel gluon exchange shown in Fig. |l[ The evolution of the vertical 
t-channel gluons to the sum over ladder diagrams is called the BFKL Pomeron which is described by the 
BFKL equation [8, 9]. 

The BFKL Pomeron splits and re-merges via the triple Pomeron vertex which was first calculated by 
Korchemsky in ref. [10] and by Bialas, Navelet and Peschanski (BNP) in ref. [11] . The large contribution 
of the triple Pomeron vertex means that Pomeron loop diagrams yield a significant contribution to the 
scattering amplitude, which is comparable to the basic amplitude of Fig. |l| (see for example estimates in refs. 
[12-14]). Hence it turns out that an accurate estimate for the scattering amplitude should take into account 
the contribution from Pomeron loop diagrams. To approach this problem requires developing an approach 
for summing over Pomeron loop diagrams. Fortunately, we can treat this theoretically in perturbative 
QCD. The sum over Pomeron loop diagrams provides the complete set of shadowing corrections to the 
basic diagram of Fig. [I] arising from hard re-scattering, which spoil the large rapidity gaps. Thus the 
formalism for summing over Pomeron loops provides the framework for estimating the contribution of 
short distance interactions to the survival probability. 




Figure 3: The special class of symmetric Pomeron loop diagrams taken into account in the sum over Pomeron loops. 

In this letter, the special class of Pomeron loops shown in Fig. [3| have been taken into account in the 
summation over Pomeron loops. Diagrams of this type can be calculated from the observation that when 
each branch of the loop in Fig. || (a) gives birth to a loop, one obtains Fig. || (b). In this way a second 
generation of loops has been introduced, and Fig. || (b) is called the N = 2 generation diagram. Fig. || 
(a) with one loop is called the N = 1 generation diagram. Likewise, a third generation of loops has been 
introduced in Fig. |3| (c) when each of the 2 simple loops at the center of Fig. |3] (b) give birth to two loops 
in the same way, so Fig. Q (c) is called the N = 3 generation diagram. Continuing in this way, the full 
spectrum of symmetric Pomeron loop diagrams is generated for all N generation diagrams. In this letter, 
an iterative algorithm is described for calculating such diagrams. 
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The solution to the BFKL equation provides the trajectory to (n, u) (n £ Z; v £ M) for the BFKL 
Pomeron as [8, 9]; 



, , 1+n \\ _ a s N c 
u(n,i>) = a s [ip(l)-Meip — h ; a s = (1) 



7T 



n represents the energy levels of the BFKL Pomeron, and v is a continuous variable which one integrates 
over when calculating Feynman diagrams. The BFKL eigenfunction falls sharply with increasing n and is 
only positive at high energy when n = 0. Hence throughout this letter which is focussed on high energy 
scattering, n = is assumed and the argument n is suppressed. Hence the BFKL Pomeron trajectory 
which is the sum over ladder diagrams of the type shown in Fig. |l] is described by the regge behavior 
s ui(u) = e u)(v)y ^ 'p^g sca ttering amplitude of Fig. || is well known and has been calculated in refs. [12-18]; 



IT 2 f°° 

A {0) (y,Sy H \Fig. § = ^-A H / dvh (v) A (u) e^ u) ^ Sm) E U E'_ U (2) 

2 11 / rn \V2+w/ r * \ 1/2-11/ 

KJ ^ K) 16(l/4 + ^) 2 ' \r 10 r 2 o) \r* 10 r%J K > 

h{v) is the integration measure which preserves conformal invariance [19, 20], A (v) is the Pomeron 
propagator in the conformal basis [19, 20] and E u is the coupling of the BFKL Pomeron to the QCD 
color dipole [18-20], in the dipole approach to proton proton scattering. Here r\2 = r\ — r 2 is the 
transverse size of the dipole and rio = r\ — tq where ro is the center of mass coordinate of the dipole. 
Syn = In (M h /Asq) (sq = lGeV/c 2 ) is the rapidity gap occupied by the heavy Higgs boson, and Ah = 
(N% — l) a s 2 1 / 4 G 1 J 2 /Sir; (where Gf is the Fermi coupling) is the contribution to the scattering amplitude 
of the process Pomeron+Pomeron — > Higgs derived in refs. [21-27], which as shown in Fig. [I] proceeds 
mainly via the intermediate top quark triangle. The observation that the BFKL eigenfunction Eq. ([!]) has 
a saddle point v = means that one can expand the exponential in Eq. (||) as 



uj (v) = uj (0) - -z/V (0) + O (u 2 ) ; w(0) = 4a fl ln2; uj" (0) = 28a s C (3) . (4) 

where £ (3) = 1.202 is the Riemann zeta function. Using this expansion the integration in Eq. (j2|) is 
evaluated by the steepest descent method which yields the result [14]; 



, , x it- Sh ^(2<r) 1/2 A„ e -«»<!<- J »»> 
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Eq. (||) is the bare scattering amplitude for the desired process of Fig. |T[ The first order shadowing 
correction is the Pomeron loop shown in Fig. ||] (a). Using the same conventions, the scattering amplitude 
is given by the expression [14]; 



A {1) (y, 5y H \Fig. | (a)) = -f A H / dvh {v) A 2 (u) e^^d {1) (u\y, 5y H ) E V E'_ V (6a) 

* J — oo 

dvih{v\)\{vi) \ dv 2 h (^2) A (z^ 2 ) |T {v\v\, v 2 ) | 2 (6b) 
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Figure 4: The triple Pomeron vertex, where lines 
indicate "reggiezed" gluons. (a) is the planar di- 
agram, and (b) is the non planar diagram. 



where y\2 = y\ — 2/2 is the rapidity gap which the loop 
fills (see Fig. |j| (a)). |r (f|fi, 1^) | 2 is the triple Pomeron 
vertex (TPV) for the splitting of a Pomeron with conformal 
variable v into two Pomerons with conformal variables v\ 
and U2, shown in Fig. [|. The explicit expression is very 
complicated and its full form and the derivation can be 
found in refs. [10, 11]. The TPV is the sum of the planar 
and non planar diagrams shown in Fig. || (a) and Fig. |4| (b), 
namely; 



1 



2tt 



T{v,v u v 2 ) = 16a 2 (- + z; 2 J ( p, (v, v u y 2 ) +^k (u, vy, v 2 ) ( X (u) - X fa) ~ X fa)) ) i (7) 



Fig. g (a) 



Fig. W (b) 



1 



In ref. [14] it was shown that the two contributions to Eq. (||) stem from the regions I; \v,iv\,iv2\ = 
{0,1/2,1/2} and region II; {iv, ^1,^2} = {1/2,0,0}, for which the TPV takes the following asymptotic 
form; 



2' 2y 2 V ' s \ N*J V4 J " W (1/2 - iv x f (1/2 - iv 2 f " ' 

1 \ l9 (4vr) 6 at 1 , N 



Eq. (|8a| ) is dominated by the planar diagram which contains singularities in region I, and the non 
planar diagram which is non singular has been thrown away. Eq. (|8b| ) is purely the contribution from the 
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non planar diagram which contains singularities in region II, whereas the planar diagram is non singular 
in region II and is not included in Eq. (|8b|). Previous publications assumed that N c — > oo, so since the non 
planar diagram has a pre- factor of 1/Nj?, it was neglected. In this letter the non planar diagram has been 
taken into account and moreover it is the dominant contribution to the TPV for region II. The non planar 
TPV leads to the contribution to the loop amplitude of Fig. [3| (a) derived in Eq. ( 10b| ) below, which is the 
dominant contribution. Hence in this letter, the remarkable property of the Pomeron loop amplitude has 
been found, that the major contribution to the Pomeron loop amplitude stems from the non planar TPV. 

Eq. (^a|) and Eq. ( pb| ) lead to two contributions to the Pomeron loop amplitude for regions I and II. 
For region I the v\,V2 integrals in Eq. ( pb| ) are evaluated by closing the contour on the upper half plane 
and summing over the residues at iv\,ivi = 1/2, and for region II the steepest descent method is used (see 
ref. [14] where the details of this calculation are explained). The loop amplitude d^ {v\y) is the sum of 



both contributions given in Eq. (9b) and Eq. (|9q), namely [14]; 
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(9a) 
(9b) 

(9c) 



Plugging Eq. (||) into Eq. (|6a| ) leads to two contributions for the scattering amplitude of Fig. || (a). 
The first part stems from Eq. ( |9a|) for region I, and the v integral is evaluated by the method of steepest 



descents. The second piece stems from Eq. (9b) for region II, for which the v integral is solved by closing the 
contour over the upper half plane and summing over the residues at iv = 1/2, and afterwards integrating 
over the rapidity variables to yield finally the two expressions [14]; 



(2ir) 1/2 a 2 dAu Mo)(y-Sy H ) / 1 „ \ 

^,fe,*»|Fig.!(a)) = 1^ - - - _ sm)f2 x (0) ( M (0) + (0) (y - 6y H) ) (10a) 

+ 2 W fe-^«)(-^j rV2(i ,_, M p < 10b » 

Eq. ( [L0a|) is the part of the scattering amplitude which leads to the renormalization of the Pomeron 
intercept u(0). Eq. fllObj ) is the dominant contribution, and is equivalent to 2 non interacting Pomerons, 
with renormalized Pomeron vertices. This can be seen from Fig. |3| (a), whereby at high energy taking out 
both branches of the loop, one is left with just two non interacting Pomeron exchanges. The same is true 
of the multiple loop diagrams of Fig. || (b) and (c), and all higher order symmetric loop diagrams, as will 
now be shown. 
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As explained above, Fig. ||| (b) stems from Fig. |3| (a) when each branch of the loop gives birth to a 
secondary loop leading to the two "second generation" of loops in Fig. |3] (b). In the same way when the 
second generation loops in Fig. |3| (b) each give birth to two loops, this leads to 4 "third generation" of loops 
in Fig. [3| (c). Continuing with this evolution, the entire spectrum of symmetric N generation diagrams 
can be generated for all N . The scattering amplitude with N generations of loops shown in Fig. || is the 
generalization of Eq. (pal) , namely [14] ; 



A (N) (y,5y H \Fig.§) = -f / dvh{u) \ 2 {v) d {N) (u\y,5y H ) E U E'_ U A H (11) 









Figure 5: The N generation diagram which stems from the simple loop giving birth to two sets of N — 1 generations 
of loops, is equivalent to the diagram of 2 N non interacting Pomerons with renormalized Pomeron vertices. 

where drpfj (v\y, Syjj) is the contribution of the iV generations of loops in Fig. ||. The way to calculate 
dm) iy\yi^VH) is based on the observation that Fig. ||| is equivalent to the simple loop diagram of Fig. || 
(a) when each branch in the loop gives birth to a set of N — 1 generations of loops. This means that to 
write the expression for d^) { u \lJi ^Uh), all that is needed is to modify the propagators for the branches of 



the loop in the expression of Eq. (6b) as; 



\{v x ) -» A(^i)d (JV -i) {vi\yi2,8yii) A(i*l) ; X(u 2 ) -> A (i/ 2 ) d(jv-i) 1 2/12) A {v 2 ) (12) 

where in the notation of Fig. ^, d(jv-i) (^1 1 2/12 5 ^Vh) labels the set of iV — 1 generations of loops on the 
right in Fig. ||| which includes the production of the Higgs, and d(jv-i) (^2 1 2/12 ) = rf(jv-i) (^lyn^VH = 0) 
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labels the set of loops on the left in Fig. |5| where there is no Higgs production. After implementing Eq. (|l2|) 
in Eq. (^) , one arrives at the following amplitude for the set of A generations of loops [14] ; 



d( N ){v\y,5y H ) = / dv\h{v\)\{ux) \ dv 2 h (v 2 ) A (u 2 ) \V {v\v\, u 2 ) 



y ry-Syn 



(13) 



x / dyi dyteW^M^Wyv-^toxd^ fa|j/i 2j 6y H ) d (N _ 1} (v 2 \y 12 ) 

J8y H JO 



Eq. (13) forms an iterative expression. Using the technique of proof by induction, the following two 
formulae for the two contributions to Eq. (]i~3| ) can be proved (see ref. [14]) which stem from the two 



asymptotic expressions for the TPV of Eq. ( pa|) and Eq. (8b); 
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Finally after inserting Eq. (fL4l) into Eq. (0), one finds the following scattering amplitude for the 
diagram of Fig. || with A generations of loops [14] ; 
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4*> (,^|F lg . 1) = LJ_^ I ^(0) ( ^ (on ^^ ))3/2 d5a) 
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A (N) (y, Sy H \Fig. g) = 29jy2 ^ y (y - 5y H ) (15b) 

j\3 p 2w(0)(3/-5j/h/2) 1 |>_1 ^\3 p 2w(0)yV 

ViV > 1. 




y 3/2 ( y _ «5y H ) 3 / 2 j [V " s d y y y 3 j 

Eq. ( [L5a| ) is the contribution which leads to the renormalization of the Pomeron intercept. Eq. ( |15b[ ) is 
the piece which is equivalent to 2^ non interacting Pomerons, with renormalized Pomeron vertices shown 
pictorially in Fig. [|. This follows from the observation that Eq. ( |15b| ) can be recast in the form; 



Af N) (y, 5y H \Fig. |) = k (n) e 2 "^ (16) 

Comparison of the RHS of Eq. (|i~6|)with Eq. (||) which is proportional to e u ^ y , shows that 
(N) (y^VHlFig. H|) is equivalent to the amplitude of 2^ non interacting Pomerons, with a pre-factor k (n) 
which contains the set of renormalized Pomeron vertices. The complete scattering amplitude is the sum of 
Eq. ( |l5a| ) and Eq. fllgED - 

Table [l] [14] shows the results up to the N = 6 generation diagram. N = refers to the basic amplitude 
of Fig. [l] without loops derived in Eq. (||). The values in the table show that for the range of energies rele- 
vant to the LHC, the scattering amplitude falls as the number of generations of loops introduced increases. 
However from N = to N = 4, for a s = 0.2 the scattering amplitudes are the same order of magnitude. 
This demonstrates the importance of taking into account the full set of Pomeron loops for the summation 
over Pomeron loop diagrams. The difference in values for different choices of a s show that the scattering 
amplitude is very sensitive to the choice of intercept. 
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A (Ar) (y = 19, Sy H = hi {M 2 H /4s )) 
a s = 0.12 a s = 0.2 


N = 


1.57 x 10~ 8 2.17 x 10~ 7 


N = 1 


2.05 x 10~ 10 2.15 x 10~ 7 


N = 2 


2.85 x 10~ 15 1.39 x 10~ 7 


N = 3 


3.40 x 10" 24 1.28 x 10~ 7 


N = 4 


8.20 x 10~ 42 2.31 x 10~ 8 


iV = 5 


8.13 x 10~ 77 5.11 x 10~ 9 


iV = 6 


1.36 x 10" 146 4.26 x 10~ 10 



Table 1: Scattering amplitudes [14] for diffractive 
Higgs production from the multiple Pomeron loop di- 
agram with N generations of loops. The mass of the 
Higgs boson is assumed to be Mr = 100 GeV, and the 
rapidity gap y between the scattering protons is taken 
to be y — 19, based on proton proton collisions at the 
typical LHC energy y/s = 14 TeV. 



The survival probability is a quantitative measure 
of the shadowing corrections to the basic diagram of 
Fig. [l| This stems from additional parton showers 
which produce inelastic scattering that destroy the large 
rapidity gaps. The contribution of hard re-scattrering 
to the survival probability can be derived from sum- 
ming over the additional parton showers that stem from 
the Pomeron loop diagrams which have been calculated 
above. In this way, the "enhanced survival probability" 
(| S^ Yx I) i s calculated by subtracting from the basic 
diagram of Fig. |TJ, the first correction from the N = 1 
generation diagram of Fig. |3| (a), and subtracting from 
this the N = 2 generation diagram of Fig. || (b) , and 
so on, and divide by the basic amplitude of Fig. [l] itself 
to yield the correctly normalized survival probability 
which leads to the formula [14]; 



(\S; 



enh I) 



1 



EM) 



jV-1 



N=l 



A {N) (y,5y H \Eq. O+Eq. fll5b|)) 
A (0) (y,6y H \Eq. ©) 



(17) 



enhanced survival probability 



0.25 



0.(15 



Fig. ^ [14] is a plot of the energy dependence of 
the enhanced survival probability. The graph shows 
a decrease of the survival probability as the energy 
increases. For the LHC range of energies, the results 
show that the enhanced survival probability is small, 
and could even be less than 1%. Both of these ob- 
servations are in agreement with the results found in 
refs. [12, 28-30]. The decrease with energy is a natu- 
ral consequence of the way the survival probability is 
denned. The survival probability measures the proba- 
bility not to have extra parton showers other than the 
basic process of Fig. ffl. Extra parton showers leading 
to unwanted inelastic scattering inevitably increase 
with the energy of the interaction, and this decreases the probability of exclusive Higgs production with 
large rapidity gaps in tact after the scattering. 




Rapidity y 



Figure 6: The enhanced survival probability (| , |) 
plotted against the rapidity separation y of the scatter- 
ing protons. a s = 0.2. 



In summary, the dominant contribution to the simple Pomeron loop of Fig. || (a) stems from the non 
planar triple Pomeron vertex shown in Fig. ||| (b), which has not been taken into account before. The 
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formula for the multiple Pomeron loop diagram for N generations of loops, has been derived in the per- 
turbative QCD approach, whereas previous methods have relied on the mean field approximation. The 
multiple Pomeron loop amplitude has two contributions which lead to the renormalization of the Pomeron 
intercept, and the diagram of 2^ non interacting Pomerons for the iV generation diagram, with renormal- 
ized Pomeron vertices. Multiple Pomeron loop diagrams, other than the simple loop of Fig. || (a) contribute 
significantly to the summation over Pomeron loops, and need to be taken into account. Using this formula, 
the contribution of short distance interactions to the survival probability has been estimated, from the sum 
over Pomeron loop diagrams, in perturbative QCD. The survival probability decreases with energy, and 
for the LHC range of energies it is small and could even be less than 1%., in agreement with the findings of 
refs. [12, 28-30]. The smallness of the enhanced survival probability show that short distance interactions 
contribute substantially. 
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discussions on the subject. This research was supported by the Fundagdo para ciencia e a tecnologia (FCT), 
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